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We use the assumption of electric-magnetic duality to express the heavy quark 
potential in QCD in terms of a Wilson Loop VK c ff (T) determined by the dynamics 
of a dual theory which is weakly coupled at long distances. The classical approxi- 
mation gives the leading contribution to W c ff(r) and yields a velocity dependent 
heavy quark potential which for large R becomes linear in R, and which for small 
R approaches lowest order perturbative QCD. The corresponding long distance 
interaction between color magnetic monopoles is governed by a Yukawa potential. 
As a consequence the magnetic interaction between the color magnetic moments of 
the quarks is exponentially damped. The semi-classical corrections to W g(r) due 
to fluctuations of the classical flux tube should lead to an effective string theory 
free from the conformal anomaly. 



1 Introduction 

In this talk we will describe how to calculate the Wilson loop W(r) determin- 
ing the spin dependent, velocity dependent heavy quark potential V qq using 
the assumption of electric-magnetic duality; namely that the long distance 
physics of Yang Mills theory depending upon strongly coupled gauge poten- 
tials is the same as the long distance physics of a dual theory describing 
the interactions of weakly coupled dual potentials C p and monopole fields Bi . 
To calculate V qq at long distances we replace W(r) by W e ff(r) a functional 
integral over the variables of the dual theorytl. Because the long distance fluc- 
tuations of the dual variables are small we can use a semi-classical expansion 
to evaluate W ff- The .classical approximation gives the dual superconductor 
picture of confinements and the semi-classical corrections lead to an effective 
string theoryp. We first review electric-magnetic duality in electrodynamics. 



2 Electric-Magnetic Duality in Electrodynamics 

Consider a pair of particles with charges e(— e) moving along trajectories 
zi(t)(z2(t)) in a relativistic medium having dielectric constant e. The tra- 
jectories Zi(t)(z2(t)) define world lines Ti(r2) running from tj to t/(t/ to ti). 
The world lines ri(T2), along with two straight lines at fixed time connecting 
Vi to 2/2 and X\ to x% , then make up a closed contour T (See Fig.l). The 
current density (x) then has the form 
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Figure 1: The loop I\ 



f{x) = e <j> dz^S{x - z). (1) 

In the usual (electric) description this system is described by a La- 
grangian 

CaU) = -\{d a A p - d (3 A a ) 2 - j a A a . (2) 

Then 

J dxC A (j) = - J dx e{ - d » Av ~ dvA ^ -ef dz^z). (3) 
The functional integral defining W(T) in electrodynamics is 

where Cgf is a gauge fixing term. 

The spin independent electron positron potential V e + e -{R,z\,z 2 ) is ob- 
tained from the expansion of i log W(T) to second order in the velocities Z\ 
and zi by the equation: 



iiogw(r) - J*' dtv e+e -(R,k u 'z 2 ) , (5) 
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where R = Z\{t) — z%{£). To higher order in the velocities, i log W(r) cannot 
be written in the above form and the concept of a potential is not defined 
because of the occurence of radiation. Eq.(|^) does not include contributions 
of closed loops of electron positron pairs to V e + e -. 
The integral (g) is gaussian and has the value 

W(T) = e-—Jr dX Jr dX = , (6) 

where D^ v is the free photon propagator. Letting e — 1 and expanding 
i\ogW(T) to second order in the velocities gives: □ 
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Furthermore the spin dependent electron positron potential V^+^L is deter- 
mined by the expectation value ((F^Y) Maxwell of the electromagnetic field in 
the presence of the external current (B]). 

In the dual description first we write the inhomogeneous Maxwell equations 
in the form: 

-^^=J„, (8) 

where G^ v is the dual field tensor composed of the electric displacement vector 
D and the magnetic field vector H: 

Gok = Hk, Gtm = tlrnnD n - (9) 

Next attach a line L of polarization charge between the electron positron 
pair. As the charges move the line L sweeps out a surface y a (a, r) bounded 
by r (the Dirac sheet) and generates the Dirac polarization tensor G? v (x): u 



yi Qy a dy^ 
da / dr — — S(x - y(cr,r)). (10) 



The current density (g) can then be written in the form: i 
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and the solution of the inhomogeneous Maxwell equations 



G>„ — d^C u — d y C^ + G^ v , (12) 
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which defines the magnetic variables (the dual potentials C M ). 

The homogeneous Maxwell equations for E and B, written in the form 



d a (nG a p) = 0, (13) 

where /i = - is the magnetic susceptibility, become dynamical equations for 
the dual potentials, and can be obtained by varying in the Lagrangian 

C C (G%) = -^G^ , (14) 

where G^y is given by ([l2]). This Lagrangian provides the dual (magnetic) 
description of the Maxwell theory (||). In the dual description the Wilson loop 
W(r) is given by 

W(T) = J-Rhill . (15) 

J VC ^J dx[Cc(GS v =0)+C GF ] 

The functional integral (JT^) is also Gaussian and has the value (||) with 
- replaced by fi. We then have two equivalent descriptions at all distances of 
the electromagnetic interaction of two charged particles. 

Note from (J2J) and ( |l4| ) that the equations 

del 9mag 

define electric and magnetic coupling constants. If the wave number dependent 
dielectric constant e — > at long distances, then g e i — > oo and the Maxwell 
potentials are strongly coupled. By contrast, g mag — > 0, and the dual 
potentials are weakly coupled at large distances. 

3 The Heavy Quark Potential in QCD 

The heavy quark potential V q q is determined by the Wilson loop W(T) of Yang 
Mills theory: 

fVAe iSYM( - A HrP exp(-ie Ldx^AJx)) 
W & = ~ JVAe^LiJ' — ( 17 ) 

(See Fig.l) As usual A^(x) = ^XaA^x), tr means the trace over color indices, 
P prescribes the ordering of the color matrices according to the direction fixed 
on the loop and Sym(A) is the Yang-Mills action including a gauge fixing 
term. We have denoted the Yang-Mills coupling constant by e, i.e., 
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(18) 



The spin independent part V(R, Z\, z%) of V q q is obtained from (17) by the 
QCD analogue of (||): 

/•*/ . . 

ilogW(T)=J dtV(R,Z!,Z2). (19) 

The spin dependent heavy quark potential V spin is-a sum of terms depend- 
ing upon quark spin matrices, masses, and momenta: a 

V Spm — Vlg AG + V T homas + y Darwin + V$S, (20) 

where the notation indicates the physical significance of the individual terms 
(MAG denotes magnetic). Each term in (20) can be obtained from a corre- 
sponding term in V by making the replacement 

((Ffiv ( z l)))Maxwell > {{Fnu {z\)))y M , (21) 

where 



Wr,{zD))YM- fDAe*s™WlrPexp[-ief r dx>'A ll (x)] ' [ ' 



and 



Fpv = d^Au - d v A M - ie[A u , A,}, (23) 

i.e. ((i r /J i/(x)))y m is the expectation value of the Yang-Mills field tensor in the 
presence of a quark and anti-quark moving along classical trajectories z\(t) 
and z*2 (t) respectively. 

The calculation of the heavy quark potential is then reduced to the evalua- 
tion of functional integrals of Yang Mills theory. Because of the strong coupling 
at long distances all field configurations can give important contributions to 
( |l7j ) and ( p2|) for large loops T and there is no simple description in terms of 
Yang Mills potentials. 
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4 The Dual Description of Long Distance Yang-Mills Theory 

The dual theory described here is a concrete realization of the Mandelstam 
't HooftH dual superconductor picture of confinement. A dual Meissner effect 
prevents the electric color flux from spreading out as the distance R between 
the quark anti-quark pair increases. As a result a linear potential develops 
which confines the quarks in hadrons. Such a dual picture is suggested bx 
the solution of a truncated set of Dyson equations of Yang Mills theory Q 
which gives an effective dielectric constant e(q) — > q 2 /M 2 as q 2 — > (M is an 
undetermined mass scale). As a consequence fi = | — > as q 2 — > so that 
the dual gluon becomes massive as is characteristic of dual superconductivity. 
However, such a truncation cannot be justified in the strongly coupled domain 
and duality in Yang Mills theory remains an hypothesis. 

On the other hand, there has been a recent revival of interest in electric- 
magnetic duality due to the work of Seiberg and WittenQ on supersymmetric 
N = 2 Yang Mills theory and Seiberg II on N — 1 supersymmetric QCD. The 
long distance physics of these models, which are asymptotically free, is de- 
scribed by weakly coupled dual gauge theories. These examples of non-Abelian 
gauge theories for which duality can can be inferred provide new motivation 
for the duality hypothesis for Yang Mills theory. 

The dual theory is described by an effective Lagrangian density £ e ff in 
which the fundamental variables are an octet of dual potentials C M coupled 
minimally to three octets of scalar Higgs fields carrying magnetic color 
charge uB. (The gauge coupling constant of the dual theory g = — ). The 
Higgs potential has a minimum at non-zero values B i which have the color 
structure 

Boi = B \ 7 , B 2 = B (— A 5 ), B 3 = B Q \ 2 - (24) 

The three matrices A7,— A5 and A2 transform as a j = 1 irreducible repre- 
sentation of an SU(2) subgroup of SU(3) and as there is no SU(3) transfor- 
mation which leaves all three Boi invariant the dual SU(3) gauge symmetry 
is completely broken and the eight Goldstone bosons become the longitudinal 
components of the now massive C M . 

The basic manifestation of the dual superconducting properties of £ e g- is 
that it generates classical equations of motion having solutions Ej carrying a 
unit of Z% flux confined in a narrow tube along the z axis (corresponding to 
having quark sources at z = ±00). (These solutions are dual to Abrikosov- 
Nielsen-Olesen magnetic vortex solutions til in a superconductor). Before writ- 
ing £ e ff we briefly describe these classical solutions. The monopole fields Bi 
have the form : El 
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B 1 = Bi (x)Xr + Bi (x) (—Xe) , 

B 2 = B 2 {x){-\ 5 ) + B 2 (x)\4, (25) 
B 3 - B 3 {x)X 2 +B 3 {x)(-X 1 ). 



We denote 



(j>i(x) = Bi(x) - iBi(x) , (26) 

and look for solutions where the dual potential is proportional to the hyper- 
charge matrix Y = 

C P = C„Y, (27) 

and where 

faix) = <t>2(x) = <j>(x), 3 (x) = B 3 (x). (28) 

At large distances from the center of the flux tube in cylindrical coordinates 
p, 8, z the boundary conditions are: 

C->- — , <j)^B e w , B 3 ^B , as p -> oo. (29) 

The non-vanishing of £>o produces a color monopole current confining the 
electric color flux. The line integral of the dual potential around a large loop 
surrounding the z axis measures this flux, and the boundary condition ( |29|) for 
C gives 

e -^l„ P a -^ =e 2-y =e Mi), (so) 

which manifests the unit of Z 3 flux in the_tube. The energy per unit length in 
this flux tube gives the string tension a: E3 

cr - 2ABl. (31) 

The field 4>(x) vanishes at the center of the flux tube. By contrast B 3 (x) 
does not couple to quarks and remains close to its vacuum value for all x. For 
simplicity in the rest of this talk we set B 3 (x) = Bq, in which case £ e ff reduces 
to the Abelian Higgs model. 

To couple C M to a qq pair separated by a finite distance we represent quark 
sources by a Dirac string tensor Gf„. We choose the dual potential to have 



the same color structure (27) as the flux tube solution. Then Gfj„ must also 



be proportional to the hypercharge matrix, 
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— YGp V , (32) 

where G^ v is given by (|To|), so that one unit of Z% flux flows along the Dirac 
string connecting the quark and anti-quark. With the ansaetze ( [32] ) and 
(p5|)- ( p8[) along with the simplification S3 (a;) = Bq , the Lagrangian density 
£ c ff (G~^J coupling dual potentials to classical quark sources moving along 
trajectories z\(t) and z%{t) assumes the form: 



CJtte _ + _ ™ A W _ sj). , (33) 



where 



and 



G^ — d^Cy — d^C^ + G?„ , (34) 



2vr , N 

.9 = — • 35 

e 

The first term in C c g is the coupling of dual potentials to quarks, the 
second is the coupling of the dual potentials to monopole fields <j>, while the 
third term is the quartic self coupling of the monopole fields. The numerical 



factors in (33) arise from inserting the color structures fl2q)- ( Bq ) in the original 
non-Abelian form of C c g. By a suitable redefinition of <fi and A the last two 
terms can be written in the standard form of the Abelian Higgs model, while 



the color factor | in the first term is a consequence of ( p7| ) and (32), which 
combined with the boundary condition ( p9| ) provides the unit of Z3 flux. 

We find from ([33|) the following values of the dual gluon mass M and the 
monopole mass : 

M 2 = 6g 2 B> , M| = ^B 2 . (36) 

The quantity g 2 /A plays the role of a Landau-Ginzburg parameter. Its value 
can be estimated by relating the difference between the energy density at a 
large distance from, the flux tube and the energy density at its center to the 

gluon condensatelij This procedure gives g 2 /X — 5. There remain two free 

2 

parameters in C e g, which we take to be a s — — and the string tension 
a. 

We denote by W e g(T) the Wilson loop of the dual theory, i.e., 
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^cff(r) = y . . (37) 

fVCpVfa'J MC o g(G^=0)+C BF ] 

The functional integral W c ff(r) determines in the effective dual theory the 
same physical quantity as W(r) in Yang-Mills theory, namely the action for a 
quark anti-quark pair moving along classical trajectories. The coupling of dual 
potentials to Dirac strings in C c g (G^„) plays the role in eq.(|3^) for W e ff(r) of 

the Wilson loop P( ,-^§ T d ^ A ^) m eq .@ for W(T). 

The assumption that the dual theory describes the long distance qq inter- 
action in Yang-Mills theory then takes the form: 

W(T) = Woff(r), for large loops T. (38) 

Large loops mean that the size R of the loop is large compared to the inverse 
of M and M$. Since the dual theory is weakly coupled at large distances 
we can evaluate W e ff(r) via a semi-classical expansion to which the classical 
configuration of dual potentials and monopoles gives the leading contribution. 
Furthermore using (|38|), we can relate the expectation value ( p2] ) of the Yang 
Mills Field tensor at the position of a quark to the corresponding expectation 
value of the dual field tensor in the effective theory: El 

«^i)»ym = l((G„M)))oS, (39) 

where 

n....(x) = 

2 



' ->Ac 



G^ix) = —£fiu\a G (x) , (40) 



and 

To obtain the spin independent heavy quark potential V(R, zi, zq) in the 
dual theory we replace W(T) by W e s(T) in eq. (|l9|) . This expresses the spin 
independent heavy quark potential in terms of the zero order and quadratic 
terms in the expansion of i log W c ff(r) for small velocities z\ and z*2- The cor- 
responding spin dependent potential in the dual theory is obtained by making 
the replacement 



{{F^( Zi ))) Y m — Tj«G^i)))cff, (42) 



4 
3 

in the expressions in eq.(f20|) for V s P in . 
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5 The Classical Approximation to the Dual Theory 

In the classical approximation all quantities are replaced by their classical 
values 



{{G lxv {x))) eS = G llv {x), ilogWgff 



dxC eS {G s 



(43) 



where G M „ and C e g (G^ v ) are evaluated at the solution of the classical equa- 
tions of motion: 



d a (d a c - d f3 c a ) = -d a G% + j$ 



■MON 



(44) 



(^-^) 2 = -^(M 2 -B o 2 ), 
where the monopole current jff ON is 

jMON = _ Ug[(j) * _ igC ^ ^ _ ( ^ + ig(J ^ , 

The boundary conditions on <fi are: 



(45) 



(46) 



ct>{x) -> 0, 



as x 



J/ (ct, t) : 



</>(x) -> B ■ 



as x 



(47) 



The vanishing of 4>(x) on the Dirac sheet y^icr, r) produces a flux tube with en- 
ergy concentrated in the neighborhood of the string connecting the quark anti- 
quark pair. Using the minimum energy solution corresponding to a straight 
line string , we evaluate i log W e g to second order in the velocities z\ , and z-i 
and obtain the spin independent heavy quark potential. 

At large separations V(R, z\, zi) is linear in R since the monopole current 
screens the color field of the quarks so that a color electric Abrikosov-Nielsen- 
Olesen vortex forms between the moving qq pair. For the case of circular 
motion, (zi ■ R — 0, z 2 — —z{), we find: 



V ^aR 



I- A 



(zi x Rf 
R 2 



as R —> oo , 



(48) 



where 



A ~ .21cr. 



(49) 



The constant A determines the long distance moment of inertia I(R) of the 
rotating flux tube: 
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lim I(R) = \{AR)R\ (50) 

At small separations the color field generated by the quarks expels the 
monopole condensate from the region between them and as R — > 0, V ap- 
proaches the one gluon exchange result, |Vd. See eq.(0). 

As the simplest application of this potential, we add relativistic kinetic 
energy terms to obtain a classical Lagrangian, and calculate classically the 
energy and angular momentum of qq circular orbits, which are those which 
have the largest angular momentum J for a given energy. We findE3 a Regge 
trajectory J as a function of E 2 which for large E 2 becomes linear with slope 
a 1 = J/E 2 = 1/8(7 (1 - A/a). Then @ gives a 1 « l/6.3cr, which is close to 
the string model relation a' — • This comparison shows how at the classical 
level a string model emerges when the velocity dependence of the qq potential 
is included. 

To calculate the spin dependent heavy quark potential we use and 
( ff3| ) to evaluate V spm (^0|) in the classical approximation to the dual theory. 
The resulting expressions are given in reference 1. Here we discuss only the 
result for the spin-spin interaction V^g ln between the color magnetic moments 
of the quark anti-quark pair. This magnetic dipole interaction is determined 
by the gradient of the Greens function G (x, x') describing the interaction of 
monopoles: 



Vss n = t^—i • - &) - (Si ■ V)(5 2 • V')G(f, x 1 ) 

6 m , i 'm 2 I I 

(51) 



G satisfies the following equation obtained from eq.(E4J) for C° 



[- v 2 +6 5 2 2 (f )]G (x , x" ) =6(2- x") , (52) 

where 4>{x) is the static monopole field. (4>(x) is real so that the monopole 
charge density j°(x) — 6g 2 (j> 2 (x)C° .) Since 4>(x) approaches its vacuum value 
Bq as x — + oo, G vanishes exponentially at large distances: 

-M\S-x'\ 

G(x,x>)^ (53) 

x^oo VK\X — X 

See eq.(||). 

The dual Higgs mechanism then produces the long distance Yukawa po- 
tential ( |53"|) between monopoles along with the linear potential (48) between 



quarks. The resulting suppression of the color magnetic interaction between 
quarks is an unambiguous prediction of electric-magnetic duality. 
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6 Fluctuations of the Flux Tube and Effective String Theory 

To evaluate the contributions to Wcff arising from fluctuations of the shape and 
length of the flux tube we must integrate over field configurations generated 
by all strings connecting the qq pair. This amounts to doing a functional 
integral over all polarization tensors G^ iv (x). Similar integrals have recently 
been carried out by Akhmedov et al. u in the case A — > oo. By changing from 
field variables to string variables ,the functional integral over G^ v (x) is replaced 
by a functional integral over corresponding. world sheets j/ M (<r, r) , multipled by 
an appropriate Jacobianpand there resultsS an effective string theory free from 
the conformal anomaly Ej Such techniques if extended to finite A could be 
applied to W e g to obtain a corresponding effective string theory . The leading 
long distance contribution to the static potential due to fluctuations of the 
string which is independent of the details of the string theory would then have 
the universal value — . EJ 

7 Conclusion 

We have obtained an expression for the heavy quark potential V q q in terms of 
an effective Wilson loop Woff(r) determined by the dynamics of a dual theory 
which is weakly coupled at long distances. The classical approximation gives 
the leading long distance contribution to W / C ff(r) and yields a velocity depen- 
dent spin dependent heavy quark potential which for large R becomes linear 
in R and which for small R approaches lowest order perturbative QCD. The 
dual theory cannot describe QCD at shorter distances, where radiative correc- 
tions giving rise to asymptotic freedom become important. At such distances 
the dual potentials are strongly coupled and the dual description is no longer 
appropriate. 

As a final remark we note that the dual theory is an SU (3) gauge theory, 
like the original Yang-Mills gauge theory. However, the coupling to quarks 
selected out only Abelian configurations of the dual potential. Therefore, our 
results for the qq interaction do not depend upon the details of the dual gauge 
group and should be regarded more as consequences of the general dual super- 
conductor picture rather than of our particular realization of it. 
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